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f— I ' Abstract 
00 

An atomic random complex measure denned on the unit disk with Normally 
distributed moments is considered. An approximation to the distribution of the 
zeros of its Cauchy transform is computed. Implications of this result for solving 
several moments problems are discussed. 
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Introduction 



Let us consider the complex random measure 

S(z) = ]T Cj S(z - z e D, Prob[i 3 e D) = 1 

where D denotes the unit disk and bold characters denote random quantities, and assume that all 
finite sets of its complex moments 

d fc = / z k dS{z) = J2 c^J, k E IN. (1) 

d i =1 

have a joint complex Gaussian distribution 

1 1 n — 1 I I 2 

g(x) = T^-^e-^^ol^-wl 
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where fj,k = E[d k ] <E(F. The Cauchy transform of S is defined as 

D ^ J= l ^ 
We are looking for the distribution of the zeros of f(z). 

Some motivations are provided in the following. We notice that, in the deterministic case, f(z) is 
equal to the derivative of the logarithmic potential 

l(*)= log no* 

Therefore the zeros of the Cauchy transform are the stationary points of L(z) i.e. the location of the 
equilibrium points in a field of force due to complex masses Cj at the points £j acting according to the 
inverse distance law in the plane. These locations were described in [7, Th.(8,2)] as a generalization 
of Lucas' theorem [7, Th.(6,l)]. 

It turns out (see e.g. [4,3] that many difficult inverse problems in a stochastic framework, can be 
reduced to the estimation of a measure S(z) from its moments d^. This is equivalent to make inference 
on p, Cj, £■ from d&. As S(z) depends in a highly non linear way on p and £ •, these are the most 
critical quantities to estimate. In [4] an approach to cope with this problem was proposed which is 
based on the estimation of the condensed density of the £j, i.e. the poles of f(z), which is defined as 

fx*-*,-) 

or, equivalently, for all Borel sets A C W 

f h(z)dz = -Y^Prob{£ 3 e A). 

A P ^ 

It can be proved that (see e.g. [5]) 

h(z) = ^Au(z) 

47T 

where A denotes the Laplacian operator with respect to x, y if z = x + iy and 

u(z) = lE\log(\fl(z-^)\ 2 ) 

P { 3=1 

is the corresponding logarithmic potential. 

However the zeros of f(^), after the logarithmic potential interpretation, convey information about 
the In the simplest deterministic case of positive Cj and p — 2, the only zero of f(z) is the 
barycenter of £1,^2 to which the weights c ^ c ^ , ei ^ C2 are associated. In the general case a qualitative 
picture of the zeros location w.r. to the £j is the following: to each couple of ^ it corresponds a 



h(z) = -E 
p 



2 



zero in a strip connecting £h,€k- The strip shape and the position of the zero in the strip is determined 
mainly by c^, c&: the smaller |c^| is w.r. to \ck\ the closest the zero is to Therefore knowledge of 
the zeros provide constraints to £j and Cj useful for estimating them (see fig. 1) . 

In the following a method to estimate the condensed density of the zeros of f(z) starting from simple 
statistics which can be computed from a sample of the moments d& is given. The idea is to interpret 
the zeros of f(z) as the poles of the Cauchy transform of a random measure S(z) associated to S(z) 
through a deterministic one-to-one transformation and to apply the same approximation technique 
used for estimating the condensed density h(z) of the poles £j in [2,1]. 



1 Pade' approximants 



In [2], in order to compute the condensed density of the poles of f(z) assuming that the joint 
distribution of a finite number of moments d , . . . , d 2p _i is Gaussian, a pencil of square random 
Hankel matrices 

U( z ) = Ui - zU 

was built from d , . . . , d 2p _i such that £ ■ are its generalized eigenvalues: 

d x d 2 ... d p 
d 2 d 3 ... dp + i 

dp dp + i . . . d 2 p_i 

This property follows from equation (1) (see e.g. [8, Sec. 7.2], [4]). But then 



U r , 



d di dp ! 
di d 2 dp 



dp i dp ... d 



2p-2 



(*) = \e |log(| f[(z - ^.)| 2 )| = E {log(| detfUo^U! - zl p }\*)} 



u 

V 

= E{log(|det[UW]| 2 |det[U 1 ]| 2 )}. 

By noticing that de^U^ 1 ]! 2 does not depend on z, and by taking the QR factorization of U(z) given 
by Gram-Schmidt algorithm we got 

h(z) = j-Au(z) = -^-ApE{\og(R kk (z)r}. 

The distribution of R 2 ^, k — 1, . . . ,p was then expanded in an uniformly convergent series of general- 
ized Laguerre functions and an explicit expression for E[log(R,kk) 2 ] and then for h(z) was found. By 
truncating the series a computable approximation of h(z) was derived as a function of the moments 
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In order to apply the methods developed in [2] we look for a new moment sequence d k 

p-i 

3=1 



and the associated measure 



= E*A* - b) 

3=1 



such that Cj are the zeros of f(z). The following theorem holds: 

Theorem 1 Let us define the first k elements d of the sequence d k for all k by 

d = T- 1 (d)e 1 (2) 

where 



T = T(d) 



Then the zeros off(z) are the generalized eigenvalues of the {p — 1) x {p — 1) pencil XJ(z) = Ui — zXJo 
where 

d 3 d 4 ... d p+ i 
d 4 ds ... d p+ 2 



do 





.. 




d 




1 


di 


d . 


.. 


, d = 


di 


, ei = 







dfc-2 • 


..d 




djk-i 








u = 



d 2 d 3 ... d T 



d 3 d 4 



dp d p+1 . . . d 2p _2 



, U 1 = 



dp+i d p+2 . . . d 2p -i 



Moreover 



p-i 

d* = E e iC*> kelN. 

.7 = 1 



(3) 



proof. Let us consider the formal random power series 

oo p oo 

F(z) = £ d k z~ k = £ c,- = 1*1 > 1 

fc=0 j=l k=0 

which can be extended to D by analytic continuation. The Cauchy transform f(z) can then be 
seen as the Pade' approximant of F(z) of orders p — l,p, denoted by [p — l,p] F (z), times z^ 1 . If 
[p — l,p] F (z) = Q p~(^ where P p (z),Q p -i(z) are polynomial of degree p,p — 1 respectively, the 
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polynomials Pfc(-z) are orthogonal w.r. to the sequence {d ,di,...} [6, Sec. 2.1] and therefore the 
zeros of Pk(z) are the generalized eigenvalues of the pencil U(z) — Ui — z\J . 

Let us define the sequence by the identity 



F^-^^d,^ (4) 
which implies that VA; G IV 

Td = ei 

The random triangular matrix T is a.s. invertible as do is a.s. different from zero. Hence d& is well 
defined. From [6, Cor. 2.7] it follows that the polynomials Qfc(z) are orthogonal w.r. to the sequence 
{d 2 , d 3 , . . .}. Therefore the roots of Q p -i(z), i.e. the zeros of f(z), are the generalized eigenvalues of 
fj(z). 

Moreover we notice that the Pade' approximant [p, p — IJf- 1 (z) to the reciprocal power series z¥~ l (z) 
is [p„x p] F ( z ) ■ Therefore the zeros of zf(z) — ]p— 1,p]f(z) are the poles of [p,p — l] F -i(2;) = Q Fp ^) • 
Because numerator has degree greater than denominator we can divide P p (z) by Q p _i(z) getting 

where Di(z) is a polynomial of degree one. But then 

p-i 



[p-2,p-l] P -i(z) = £ 



C 3 



and from eq.4, we get d k = Y?j=i CjCj- a 



2 Approximate distribution of the moments of the associated measure 

In order to derive the joint density of d , . . . , d n _x we need the following lemma 
Lemma 1 The transformation $ : d — > d defined by 

d = T- 1 (d)e 1 (5) 
is an involution i.e. $($(d)) = d. 

proof. The thesis follows by the Toeplitz structure of the matrix T. □ 
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Theorem 2 The joint density of do, . . . , d n _ x is given by 



g(y) = 



(7r<7 2 j/g)» 



e ^ 



where 



and 



(-fc+2) 



det 



.j_^fe(fc+l)/2jy(-fc+2) 



VlV2 ■■■ Vn-1 Vn 
y yi ... y n -2 Vn-1 

■■•2/0 yi 



proof. If n = 2p the joint density of d , . . . , d„_x is 

1 



1 ELoi^-^i 



7TO" 



2\n 



e 



By Lemma 1, if y — <&(x), x,y E(F n , the complex Jacobian of the transformation $ is 

2/o 2n 2/o 2n 

because n is even. But then, the joint density of d , . . . , d n _i is 

1 



g(y) 



-e ° 



-Efciol[* _1 (w)]*-wl 2 



(™ 2 y 2 )' 

The rest of the thesis follows by [6, pg.106]. □ 

Theorem 3 In the limit for a — > 0, g(y) tends to a Dirac measure centered in $(//). 
proof. By noticing that g(y) can be written as 

9(y) = J S(y-^(x))g(x)dx 

in the limit for a — > we get, for every continuous function if)(y) with compact support 

Jim J g(y)ip(y)dy = J™ ^ J S(y - $(x))g(x) 4>(y)dxdy = 

lim J g(x) (^J 5{y — 3>(x)) ip(y)dyj dx = lim J g(x)ip(^(x))dx = 
J S(x - fi)ij>($(x))dx = ^($(^)). □ 



As a consequence of this theorem we can approximate g(y), in the limit for o — > 0, by a Gaussian 
centered in $(//) with height 



i.e. 

// 2n u 2 

- (ira 2 ) n 



3 Approximate condensed density of the zeros 



In view of the result of the previous section, we can use the theory developed in [2,1] to find an 
approximation of the condensed density of the zeros of f(z). More precisely starting from new data 
do, ... , d n _i that are approximately distributed as a Gaussian, we can compute the QR factorization 
of the pencil U(z). The squared diagonal elements R| fc , k = 1, . . . ,p — 1 are then (conditionally) 
distributed as quadratic forms in Normal variables and then their distribution can be approximated 
by a random variable whose density fk(y) admits a uniformly convergent Laguerre series ( [2, Th.l]) 



fk(v) = h$ t)V 



+ 



E i&L m (y/T k ,a k ). 



m=l 



But then E[\og(Tlkk) 2 } can be approximated by taking the expectation of each term of the series 
which can be computed in closed form giving 

E {log \R kk (z)\ 2 } » ^(^[log /^) + *(a k (z))] + 
- ™ T(a k (z) + h) (0k(z)\ h n ofS T , . N , s. i 

E CM E v r ^ 7^ iog^(^) + + >0 

m =i ^=0 r (M*)) Vw)/ 

where c^ m are the coefficient of the Laguerre polynomials, b$(z) are the coefficient of the Laguerre 
expansion, a k (z), (3 k (z) are the parameters of the Gamma density which is the first term of the series 
for fk(y), and is the logarithmic derivative of the Gamma function. 

The estimate of the condensed density of the zeros of f(z) is then obtained by truncating the series 
after the first term, replacing a k (z) by convenient estimates, and (3 k (z) by j3 > which turns out to 
be a smoothing parameter ([2, Sec.3]): 



p-i i 
h(z) oc E A * 

k=i V 



a 2 (3 



where A is the discrete Laplacian, and R kk {z) is obtained by the QR factorization of a realization 
of the pencil \J(z). 
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4 Simulation results 

To appreciate the goodness of the approximation to the density of R| fc provided by the truncated 
Laguerre expansion, for n = 74, N = 10 6 independent realizations k — 1, . . . , n, r — 1, . . . , N 
of the r.v. d& were generated by the equation 

4 ] = s k + 4\ «* = 5># 

3=1 

where ej^ are independent realizations of zero mean Gaussian variables with standard deviation 
a = 0.2, p* = 5 and 

-0.1-i27r0.3 -0.05-i27r0.28 -0.0001+i27r0.2 -0.0001+i27r0.21 -0.3-i2tt0.35" 

c= [6,3,1,1,20]. 

The matrix T^ r ^ based on was formed and the linear system (2) was solved for the new moments 
d~k\k = l,...,n, r — 1, . . . , N. The (p — 1) x (p — 1) matrices Uq\ based on d^ were 
computed. The matrix Ui ^ — zUq ^ with z = cos(l) +i0.8 was formed, its QR decomposition and the 
first 10 empirical moments jj were computed. Estimates of the first 10 coefficients of the Laguerre 
expansion were then computed by ([9]) 

A if R ^ ~ H 

ak = ~ ro' > j k = : 

72 - 7i 7i 

&? ) H-i) h rWE(-i) J ( lj H v 70 = 1, h = i,...,io 

The one term and ten terms approximations of the density were then computed and compared with 
the empirical density of R| fc for k — 1, . . . , p. The results are given in fig. 2. In the top left part the real 
part of one realization of the transformed data dk and transformed signal Sk are plotted. In the top 
right part the L 2 norm of the difference between the empirical density of R| fc , k — 1, . . . ,p computed 
by MonteCarlo simulation and its approximation obtained by truncating the series expansion of the 
density after the first term and after the first 10 terms is given. In the bottom left part the density 
of R|fc, k = 35, approximated by the first term of its series expansion and the empirical density are 
plotted. In the bottom right part the density of R| fc , k = 35, approximated by the first 10 terms 
of its series expansion and the empirical density are plotted. It can be noticed that the first order 
approximation is quite good even if it becomes worse for large k. 

To appreciate the advantage of the closed form estimate h(z) with respect to an estimate of the 
condensed density obtained by MonteCarlo simulation an experiment was performed. iV = 100 
independent realizations of the r.v. generated above were considered. 
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An estimate of h(z) was computed on a square lattice of dimension m = 100 by 



N p 



r=l fc=l 



J2 



a 2 (3 



+ 1 



where R^ r \z) is obtained by the QR factorization of the matrix — zUq T> . In the top part of 
fig. 3 the estimate of h(z) obtained by Monte Carlo simulation is plotted. In the bottom part the 
smoothed estimates h(z) for o = 0.2 and = 5n based on a single realization was plotted. 

We notice that by the proposed method we get an improved qualitative information with respect to 
that obtained by replicated measures. This is an important feature for applications where usually 
only one data set is measured. 



r(r) 
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Fig. 2. Top left: real part of the transformed signal (solid) and transformed data (dotted) with a = 0.2; 
top right: L2 norm of the difference between the empirical density of R|fc,& = 1, . . . , 35 computed by 
MonteCarlo simulation with 10 6 samples and its approximation obtained by truncating the series expansion 
of the density after the first term (dotted) and after the first 10 terms (solid); bottom left: density of 
R| fc , k = 35, approximated by the first term of its series expansion (solid), empirical density (dotted); 
bottom right: density of R| fe , k = 35, approximated by the first 10 terms of its series expansion (solid), 
empirical density (dotted). 
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Fig. 3. Top: Monte Carlo estimate of the condensed density when a = 0.2 based on 100 samples; bottom: 
estimate of the condensed density by the closed form approximation with (3 = 14 based on one sample. 
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